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ABSTRACT 



We report the results of a high-energy multi-instrumental campaign with INTEGRAL, RXTE, and Swift of the recently discovered 
INTEGRAL source IGR J19294+ 1 8 16. The SwiftfKRT data allow us to refine the position of the source to RAj 2000 = 19' 29"' 55.9 s 
Decj2ooo=+18° 18' 38.4" (±3.5"), which in turn permits us to identify a candidate infrared counterpart. The Swift and RXTE spectra 
are well fitted with absorbed power laws with hard (T ~ 1) photon indices. During the longest Swift observation, we obtained evidence 
of absorption in true excess to the Galactic value, which may indicate some intrinsic absorption in this source. We detected a strong 
(P=40%) pulsation at 12.43781 (+0.00003) s that we interpret as the spin period of a pulsar. All these results, coupled with the 
possible 1 17 day orbital period, point to IGR J19294+ 1816 being an HMXB with a Be companion star. However, while the long-term 
INTEGRAL/IB1S/1SGR1 18-40 keV light curve shows that the source spends most of its time in an undetectable state, we detect 
occurrences of short (~ 2000 - 3000 s) and intense flares that are more typical of supergiant fast X-ray transients. We therefore cannot 
make firm conclusions on the type of system, and we discuss the possible implications of IGR 1 1 9294+ 1816 being an SFXT. 
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1. Introduction 



Accretion, accretion disks ; Stars: individual: IGR 119294+1816, IGR Jl 1215-5952, IGR 



The INTErnational Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) has permitted a large number of new sources to 
be discovered. Amongst the ~ 250 new source^ INTEGRAL 
has unveiled two new or poorly-known types of high mass X- 
ray binaries (HMXB): the very absorbed supergiant HMXBs, 
and the supergiant fast X-ray transients (SFXT). While all types 
of HMXBs (including those hosting a Be star, Be-HMXBs) are 
powered by accretion of material by a compact object, under- 
standing the nature of a system is very important for studying 
the evolutionary paths of source populations, and more globally, 

the evolution of the Galaxy in terms of its s ource content . 

IGR J 19294+ 18 16 was discovered by iTuerler etail d2009l) 
who reported activity of this source seen with INTEGRAL 
during the monitoring campaign of GRS 1915+105 (e.g. 
IRodriguez et al.Ll2008al) . Archival Swift data of a source named 
Swift J1929. 8+1818 allowed us to give a refined X-ray posi- 
tion of RAj2ooo= 19* 29'" 55.9 s and Dec J2 ooo=+18 18' 39" 
(±3. 5") which we used to locate a possible infrared counter- 
part (IRodriguez et al.L 120091) . We also suggested that the Swift 
and INTEGRAL sources are the same and that activity from 
IGR J19294+1816 had been seen with Swift in the past. The 
temporal analysis of the XRT light curve sho wed a possi- 
ble pulsation at 12.4 s (IRodriguez et al.L 120091) . Analysis of 
Swift/BAT archival data revealed that the source had been de- 
tected on previous occasions with a periodicity at 117.2 days 
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1 see http://isdc.unige.ch/~rodrigue/html/igrsources.html for an up to 
date list of all sources and their properties! 



(ICorbet & Krimml |2009l) that we interpreted as the orbital pe- 
riod of the system. 

Soon after the discovery of the source with INTEGRAL, we 
triggered our accepted RXTE and Swift programmes for follow- 
up observations of new INTEGRAL sources (PI Rodriguez). 
A preliminary analysis of the real-time RXTE data allowed 
us to confirm the puls ations in the signal from the source 
(IStrohmaver et al.L 12009) at a barycentred period of 12.44 s indi- 
cating that this object hosts an accreting X-ray pulsar. Here we 
report the detailed analysis of the INTEGRAL, Swift, and RXTE 
observations. We begin this paper by detailing the procedures 
employed for the data reduction. We then describe the results (re- 
fined position, X-ray spectral, and temporal analyses) in Sects. 
3, 4, 5, and 6, and discuss them in Sect. 7. 

2. Observations and data reduction 

The log of the Swift and RXTE observations of 
IGR J19294+1816 is given in Table Q] The LHEASOFT 
v6 . 6 . 2 suite was used to reduce the Swift and RXTE data. The 
RXTE/PCA data were reduced in a standard way, restricting it to 
the times when the elevation angle above the Earth was greater 
than 10°, and when the offset pointing was less than 0.02°. 
We accumulated spectra from the top layer of proportional 
counter unit 2 (PCU 2), the only one that was active during all 
the observations, and therefore time-filtered the data for PCU 
2 breakdown. The background was estimated using the faint 
model. In addition, we produced event files from Good Xenon 
data with make_se, from which we extracted 2~ 12 s (~ 250 fis) 
light curves between 2 and 20.2 keV (absolute channels 5 to 48 
in order to limit the instrumental background at low and high 



2 



Rodriguez et al.; IGR J19294+1816 as seen with INTEGRAL, RXTE, and Swift 



Table 1. Journal of the Swift and RXTE dedicated observations. 
Obs. S3 and X5 were simultaneous with an INTEGRAL obser- 
vation. 



Satellite 


Date start 


Exposure time 


Name 




(MJD) 


(s) 




Swift 


54443.051 


7959 


SI 


Swift 


54447.013 


3457 


S2 


RXTE 


54921.315 


2600 


XI 


RXTE 


54921.708 


3344 


X2 


RXTE 


54922.813 


3344 


X3 


RXTE 


54923.730 


1440 


X4 


Swift 


54925.814 


2565 


S3 


RXTE 


54925.826 


3360 


X5 



energies) with seextrct. Light curves with 1-s time resolution 
were obtained in the 2-60 keV range from standard 1 data. 
These light curves were then corrected for the background and 
used for the temporal analysis. 

The Swift/XRT data were first processed with the 
xrtpipeline vO.12.2 tool to obtain standardly filtered 
level 2 event files. Images, light curves, and spectra were 
extracted within XSELECT V2 . 4a. Source light curves and 
spectra were extracted from a circular region of 20 pixels radius 
centred on the most accurate source position available. We 
tested different regions of extraction for the background light 
curves and spectra, and obtained no significant differences. 
For the remainder of the analysis, the background region is 
a circular region of 60 pixels radius centred at RA=19'' 29" 1 
41.8 s and Dec=+18° 21' 11.5". The ancillary response file was 
obtained with xr tmkarf 0.5.6 after tak ing the exposure map 
into account (see lRodriguez et al 1 l2008bL for more details). 

We also analysed all INTEGRAL public pointings aimed at 
less than 10° from the position of IGR J19294+1816 plus all 
private data of our monitoring campaign on GRS 1915+105. 
The Offline Scientific Analysis Software (OSA) 
v7 . was used for the INTEGRAL data reduction. We further 
filtered the list of science windows (SCW) to those where the 
good exposure time of ISGRI was greater than 1000 s. After 
omitting the bad pointings (e.g. first or last pointings before 
the instruments are switched off due to passage through the 
radiation belts) this resulted in a total of 1476 good SCWs. 
These pointing s were blindly s earched for the presence of 
the source (see IRodriguez et al.l d2008al) for the procedure of 
data reduction and source searching). An updated catalogue 
containing IGR J19294+ 18 16 with the most precise Swift 
position was given as an input to the software. Mosaics were 
accumulated on a revolution basis, but this results in mosaics 
of very different exposure times. Note that even in revolution 
788, w here we confirm the detection reported by iTuerler et alJ 
d2009l) . the source is too dim for further spectral analysis. 



3. IGR J1 9294+1 81 6 as seen with INTEGRAL 

IGR J19294+1816 w as discovered on Ma rch 27, 2009 (revolu- 
tion 788, MJD 54917 lTuerler et alll2009h . A re-analysis of the 
consolidated data allowed us to confirm a source detection dur- 
ing this observation. The source was detected at 10.7 cr between 
18 and 40 keV, with a flux of about 17 mCratQ. A re-analysis 
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2 The count rate to Crab conversion is based on an analysis of the 
Crab observation performed during rev. 774, with 1 Crabig_40keV = 
181.1 cts/s. 



Fig.l. INTEGRAL/IBIS 18-40 keV light curve of 
IGR J19294+1816 showing the flare on MJD 53861. 



of the previous observation on March 21, 2009 (rev. 786, MJD 
5491 1), led to a detection in the 18-40 keV energy range with a 
significance of 4.5 cr and a flux of 7.2 mCrab. IGR J19294+1816 
was not detected on April 4, 2009 (rev 790, MJD 54925) with a 
18-40 keV 3-cr upper limit of about 5 mCrab. The source was 
not detected in any of the mosaics from the previous revolution. 
However, given the very different exposure times we used in 
building the mosaics and the periodicity of t he source activity 
seen with SwiftfBKT dCorbet & Krimml 120091) . this may simply 
indicate that the periods of activity of the source are short. 

Examining these data on a SCW basis shows that the source 
was detected on a few occurrences. In particular, it was clearly 
detected during an observation in rev. 435 (MJD 53861) and 
during another pointing in rev. 482 (MJD 54004). In the for- 
mer, IGR J 19294+ 18 16 was detected in SCW 9 (19.17-19.77 h 
UTC) with a 18-40 keV flux of 33.7 mCrab. It was not detected 
in the following SCW (19.85-20.50 h) with a 3-cr upper limit 
of 11.6 mCrab. No pointings are available before SCW 9. The 
18-40 keV light curve of this epoch is reported in Fig.Q] In rev. 
482, the source was detected in SCW 51 (1-05-2.03 h UTC) at a 
18-40 keV flux of 22 mCrab (5.4 cr). It was not detected in the 
previous and following SCWs with respective 3-cr upper limits 
of 11.0 and 13.9mCrab. 

4. Refining the X-ray position, and the search for 
counterparts 

We combined all Swift observations to create a mosaic image of 
the field. As shown in Fig.|2j a single bright X-ray source can be 
seen wit hin the ~ 5' 90% erro r circle of the INTEGRAL position 
given in ITuerler et al.l d2009l) . This source lies at (as obtained 
with xrtcentroid)RA J2 ooo=19 / ' 29'" 55.9 s and Dec J2 ooo=+18 
18' 38.4" with an uncertainty of 3.5". Note th at this is slightly 
differe nt from the preliminary position given in IRodriguez et al.l 
(120091) . but the two positions are clearly compatible. Given the 
addition, here, of an additional Swift pointing to obtain these co- 
ordinates, we consider this position as the most precise one. 

We searched for counterparts in catalogues available online 
such as the 2 Microns All Sky Survey point sour ce and extended 
sourc e catalogue^ (2MASS and 2MASX, ISkrutskie et all 
120061) . the National Radio Astronomy Observatory VLA Sky 
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Fig. 2. SwiftfKRT mosaic of the field around IGR J19294+1816. 
The circle represents the error box from INTEGRAL. 



Table 2. Infrared magnitudes dereddened with different values 
of the absorption. 



N H 


A v 


Jder. 


Hder 


K s-C lcr. 


J-H 


H-K s 


xlO 22 cm- 2 














2.0 


11.20 


11.32 


11.04 


10.86 


0.28 


0.18 


2.73 


15.29 


10.14 


10.32 


10.40 


-0.18 


-0.08 


4.0 


22.40 


8.08 


9.09 


9.53 


-1.01 


-0.44 



Survey (NVSS, ICondon et ail 1 19981) . and the catalogue of the 
US Naval Observatory (USNCQ). We also made use of im- 
ages from the Second Palomar Observatory Sky Survey obtained 
through the STScI Digitized Sky Survey (DSS0). Finding charts 
from 2MASS IR and DSS are shown in Fig.fJ] 

A single infrared counterpart is listed in the 2MASS cat- 
alogue: 2MASS J19295591 + 1818382 has magnitudes J=14.56 
+0.03, H=12.99 +0.02, and K s =12.11 ±0.02. No optical or ra- 
dio counterpart was found in any of the other catalogues, al- 
though a faint object might be present in the "infrared" (7000- 
9700A) POSS II image (Fig. [3j. The infrared magnitudes were 
dereddened using Aa = C(A)xAy, with C(A) the wavel ength 
dependent coefficients o b tained by ICardelli et al.l d!989l) (see 
Table 4 in iRahoui et all 2008L for the valu es adopted), and 
Ai/=5.6xl0 22 x N H dPredehl & Schmittl[T995l) . The value of the 
total absorbing column density along the line o f sight was ob- 
tained from the LAB surv ey of Galac t ic HI dKalberla et all 
120051) . and the CO survey of iDame et ail (1200 ll) for the value of 
molecular H 2 . Therefore, N H = N m + 2N Hl = 1.23 + 2 X 0.75 = 
2.73 x 10 22 cirT 2 . We also used the value of Nh obtained through 
the X-ray spectral fits to the data, Nex,bbody=2.0xl0 22 cirr 2 , and 
Nhx.po= 4.0xl0 22 cm -2 (see Sec. |5J. The dereddened magni- 
tudes are reported in Table [2] 

4 http://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/icas 

5 http://archive.stsci.edu/cgi-bin/dss_form 



5. X-ray spectral analysis 

5. 1 . Correction of the Galactic ridge emission in the PC A 
spectra 

Since IGR J19294+1816 is faint and located in the Galactic 
plane (its Galactic coordinates are / = 53.5400°, b = 0.1150°), 
the Galactic ridge can contribute a significant amount of the flux 
collected by the PCA. We estimated the lev el of emission using 
the results from IValinia & Marshalll d 1998b . and first corrected 
the spectra from the Galactic rid ge emission in a similar way 
as done for IGR J19140+0951 bv lPrat et alj d2008l) . i.e. adding 
the Galactic ridge spectrum to the instrumental background 
spectrum. 

A simultaneous fit to the S3 and X5 spectra showed a large 
offset in the normalisations of the spectra: when including a con- 
stant in the fit and freezing it to 1 for the S3 spectrum it resulted 
in a low value of ~0.45 for the X5 spectrum. This might indicate 
that either the source varied during the non-simultaneous part 
of these observations, or that we over-corrected the emission 
of the Galactic ridge in the PCA spectrum. We consider this 
last hypothesis as the most likely given that the spec tr um of 
the Galactic ridge obtained from IValinia & Marshalll dl998l) 
is an average spectrum covering the central regions of the 
Galaxy, i.e. -45° < I < 45° and -1.5° < b < 1.5°, and that 
IGR J19294+1816 is slightly outside this region. Although 
the source is located in the Galactic ridge, we can expect 
local variations in the ridge, and/or slight changes (lower 
intensity for example) at high longitude (as also mentioned in 
IValinia & Marshalj|1998l) . 

In a second run, we did not correct the X5 spectrum 
for the contribution of the Galactic ridge, and we did not 
include a normalisation constant between the S3 and X5 
spectra (we therefore made the assumption that Swift and 
RXTE are perfectly cross calibrated). The spectra were fit 
with an absorbed power law and th e best-fitting model for th e 
Galactic ridge spectrum obtained bv lValinia & Marshalll dl998l) . 
A multiplicative constant precedes this spectral component 
whose parameters are frozen. The spectral model is therefore 
phabs*powerlaw+constant* (wabs* (raymond+powerlaw)) 
in the XSPEC terminology. We initially fitted the S3 spectrum 
alone (with the constant fixed at 0) to obtain the parameters of 
the source. Then, we added the X5 spectrum and only left the 
normalisation constant to the X5 data free to vary while it was 
held at for the S3 data. Therefore, we assume that the average 
spectral parameters obtained in the -45° < / < 45° are still valid 
at the position of IGR J19294+1816 , and that only the overall 
flux of the ridge emission changes. The most precise value 
obtained for the constant is 0.85+0.04 (at 90% confidence). 
In subsequent analyses, we used this value to correct all PCA 
spectra according to: 



BGD tot = BGD instl + 0.85 x Ridge 



5.2. Results of the spectral analysis 



(1) 



The Swift/XRT spectra were rebinned so as to have a mini- 
mum of 20 cts per channel, except for the channels below about 
1 .2 keV that were grouped in 2 bins: one below 0.2 ke V that was 
not considered further for the spectral fitting, and one between 
0.2-1.2 keV to obtain an additional spectral point to better con- 
strain the absorption column density Nh- The Swift spectra were 
fitted between ~ 0.8 and ~ 8 keV, while the RXTE spectra were 
fitted between 3 and 25 keV. We fitted each individual spectrum 
alone in XSPEC 11.3. 2ag, except the spectra from observations 
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Fig.3. 1.14'x0.73' POSS II optical (top) and 2MASS infrared (bottom) images of the field around IGR J19294+1816. The black 
circle is the Swift error box. a) " Blue", b) "Red", c) "infrared", d) J filter e) H filter, f) K s filter. 



S3 and X5 that were fitted simultaneously. Different spectral 
models were used to fit the spectra. The Swift spectra are well 
fitted using either an absorbed blackbody or an absorbed power 
law. Only the latter model provides a good fit to the RXTE spec- 
trum, and joint Swift and RXTE spectra. Considering the detec- 
tion at energies greater than 20 keV with INTEGRAL and Swift 
, this suggests that the true origin of the emission is probably 
not from a blackbody. However, we report the parameters of 
the blackbody together with those obtained with the power-law 
model in Table [3] Observations S3 and X5 were simultaneous 
with an INTEGRAL observation, but the source was not detected 
by the latter satellite with an 18-40 keV 3-cr upper limit of 
5 mCrab with IBIS. 

Since the PCA is well calibrated only above 3 keV, it is not 
suitable for proper measures of Nh, especially in such a dim 
source. In fact, leaving all parameters free to vary during the 
spectral fits leads to poorly-constrained values for this parameter 
(N H < 5.4 x 10 22 cm -2 in Obs. XI). We therefore froze the value 
of Nh to the value returned by the fit to the closest Swift data 
(Obs. S3 which, simultaneously fitted with X5, leads to 3.4 x 
10 22 cm -2 for a power-law model). 

6. X-ray temporal analysis 

Although the fir st hint for the pre s ence o f a pulsation came from 
Swift/XRT data dRodriguez et all 120091) . we focus here only on 
the data from RXTE/PCA which is the instrument dedicated to 
the timing analysis of astrophysical sources. Also, we only con- 
sider here Obs. XI and X3 which both have 2 PCUs on and 
exposure times longer than 2.5 ks, in order to increase the sta- 
tistical significance of the pulse. We produced power spectral 
density (PSD) distributions between 0.00195 and 1024 Hz with 
powspec vl . (8 averaging all sub-intervals of 512 s long of each 
of the two observations. The continuum of the 0.00195-1024 Hz 
PSD is equally well represented by a model consisting of a con- 
stant (to account for the contribution of Poisson noise), and ei- 
ther a power law (of index F = -1.35 + 0.08) or a zero-centred 
Lorentzian (with a width of 1 .6 + 0.4 x 10~ 2 Hz). Fig.|4]displays 
the white-noise corrected PSD restricted to the 0.00195-1 Hz 




0.01 0.1 1 

Frequency (Hz) 



Fig. 4. 0.00195-1 Hz white-noise corrected power spectral den- 
sity distribution of RXTE Obs. 1 & 3 obtained from the 2- 
20.2 keV high-resolution light curve. 



range. Above ~ 0.5 Hz, the PSD is compatible with white noise. 
A strong peak at around 0.08 Hz and its first harmonic at about 
0.16 Hz can be seen above the PSD continuum (Fig. |4). The 
main peak (at 0.08 Hz) is contained within a single frequency 
bin which is indicative of a coherent pulsation. However, we 
point out that the bottom of this thin peak is broadened. This 
broad feature can be fitted with a Lorentzian of centroid fre- 
quency 0.078 + 0.002 Hz, and a width of 0.020 + 0.008 Hz. 
It has a (raw) power of 11.1% RMS. We also note the pres- 
ence of a possible feature at ~ 0.035 Hz, yet the inclusion of 
a second Lorentzian to the model does not improve the fit in a 
significant way. To study the properties of the coherent pulsa- 
tion, we barycentred the 1-s light curves with barycorr vl . 8 
using the orbital ephemeris of the satellite, and we corrected 
the light curves for the background. The pulsation period was 
searched for using a Lo mb-Scargle periodogram as described in 
iPress & Rvbickl (1 19891) with the errors calculated from the peri- 
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Table 3. Results of the X-ray spectral fitting. The errors are at the 90% confidence level. 
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Obs N H T kT bb Fi-iowabs. xi* 





(xlO 22 cm" 2 ) 




(keV) 


(erg cm 2 s 1 ) 


(DOF 1 ) 


SI 


4 


0.9+0.2 




3.6 xlO -11 


1.0 (66) 




2.0±0.4 




2.1±0.2 


3.0X10" 11 


0.86 (66) 


S2 


3 4 +L0 

J - -0.8 


1.0+0.3 




3.8 xlO-" 


0.78 (34) 




1.8±0.5 






3.0 xlO" 11 


0.74 (34) 


XI 


3.4 (Frozen) 


i 17+O.I4 
1,1 '-0.11 


2.65 xlO- 11 


0.94 (25) 


X2 


3.4 (Frozen) 


1 +0 3 




0.9 xlO-" 


0.49 (10) 


X3 


3.4 (Frozen) 


1.2+0.1 




2.6x10"" 


0.91 (25) 


X4 


3.4 (Frozen) 


1.1+0.4 




1.1x10"" 


0.93 (12) 


S3+X5 


3 4 +LS 


1 2 +03 




0.74x10-" 


0.86 (8) 



* Reduced^ 2 . 

1 Degrees of freedom. 



150 



100 



50 




.Ji.,i,i...i...iL.iu,,l., iJuA^j] iUktUiilu itJlLtofl iWikAojUi^W 

j i i i i i < I i i i i i i i i I 



10 

Period (s) 



100 



Fig. 5. Lomb-Scargle periodogram of RXTE Obs. 1 & 3. The 
prominent peak is found at a period of ~ 12.44 s. The insert 
shows the light curve folded at 12.44 s. 



odogram foll owinglHorne & Baliunasl (Il986l) (see also the orig- 
inal work of Kovacs , 1 19811) . This is the same method used by 
IRodriguez et al.l ( 20061) to find the pulse's most precise period 
in IGR J16320-4751. The light curves were folded with efold 
vl . 1. We combined Obs. XI & X3 to obtain the most precise 
measure of the pulsation. A very prominent peak is visible in 
the periodogram (Fig. \5}, together with its first harmonic. We 
derived a pulse period of P=12.43781(3) s (0.080400 Hz). The 
light curve, folded at a period of 12.44 s, is shown in Fig. [5] as 
an insert. The pulse fraction was calculated in a standard way 
(P = , see e.g. IRodriguez et all 120061) . Without taking 

the emission of the Galactic ridge into account, we estimate a 
pulse fraction P — 28 + 2.3%. Given the source's location in the 
Galactic plane which has a non-negligible contribution to the to- 
tal flux, this value underestimates the true pulse fraction, and 
should be considered as a lower limit. According to the results 
of the spectral analysis presented in the previous section, we ob- 
tained the "true" background following Eq. Q] which includes 



the estimated contribution from the Galactic ridge at the posi- 
tion of IGR J19294+1816. The light curves corrected for this 
total background yields P = 40 + 4.4%. 



7. Discussion and conclusions 

We have reported here the results obtained from a multi- 
instrumental campaign dedicated to the X-ray properties of new 
INTEGRAL sources. The refinement of the X-ray position pro- 
vided by the Swift/XRT observations enabled us to identify a 
possible counterpart at infrared wavelengths. The differences 
of dereddened magnitude s (Table [3 do n ot lead to any of the 
spectral types tabulated in lTokunagal(l2000h . With N H =2.0xl0 22 



cm , J-H=0.28 would indicate an F7 V or F8 I star. However 
the value of H-K s is inconsistent with both possibilities. With 
Nh=2.73x10 22 cm- 2 , J-H seems too high for any spectral type, 
although we remark a marginal compatibility (at the edge of 
the errors on the magnitude) w ith an 09.5 V star (J-H=-0.13, 
H-K g =-0.04 lTokunagall2000l) . This value is, however, a mea- 
sure of the interstellar absorption through the whole Galaxy. 
Since IGR J 19294+ 18 16 probably lies closer than the other 
end of the Galaxy, it is likely that the interstellar absorption 
along the line of sight is lower. In particular, we note that with 
N H =2.5xl0 22 cm 2 , we obtain J-H=-0.04 and H-K s =0.05, 
which is very close to t he values tabulate d for a B3 I star 
(J-H=-0.03, H-K s =0.03 iTokunagal l2000h . In that case, the 
source would lie at a distance d > 8 kpc. 

The X-ray behaviour of the source is indicative of an HMXB: 
the X-ray spectra are power law like in shape with a hard photon 
index. One of the spectra (Obs. SI with the longest exposure) 
shows evidence for absorption in clear excess to the absorption 
on the line of sight, which may indicate some intrinsic absorp- 
tion in this source. However, the evidence is marginal in the other 
spectra, which could suggest that the intrinsic absorption varies 
in this system, as has been observed in a number of HMXBs. A 
long-term periodicity was revealed in the Swift [RPR data w hich 
confirms the binarity of the source (ICorbet & Rrimmi r2009). We 
clearly detect an X-ray pulse at a period of 12.44 s, indicating the 
presence of an X-ray pulsar. We estimate a pulse fraction P - 
40%. Apart from millisecond X-ray pulsars, found in systems 
that are at the end of the evolutionary path of X-ray binaries (and 
are LMXBs), pulsars are young objects that are usually found in 
HMXBs. We also remark that the PSD of IGR J19294+1816 
shows a broadening at the bottom of the coherent pulsation. 
Sidelobes and other noise features around co herent pulsation 
signals have been seen in other HMXBs (e.g. iKommers et all 
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Fig. 6. Most recent version of the "Corbet diagram" including 
all sources detected with INTEGRAL (Bodaghee et al. 2009 
in prep.). Squares represent Be-HMXBs, circles Supergiant- 
HMXBs (Sg-HMXB), circles in squares are SFXTs, filled sym- 
bols represent sources discovered by INTEGRAL. The triangles 
are HMXBs of unknown nature (including IGR J19294+1816). 
The positions of IGR J 1 9294+ 1 8 1 6 and of the two SFXTs lying 
in the Be region of the plot are highlighted. 



[T9^lKauretail.l2007l in, respectively, 4U 1626-67, and XTE 
J01 1 1.2-7317). Such features can be produced by artificial ef- 
fects (such as the finite length of the time intervals used to make 
the PSD), or they can be real if, e.g., a q uasi-periodic osc i llation 
(QPO) beats wit h the coherent signal dKommers et all 119981 : 
lKauretalll2007h . In the case of IGR J19294+1816, a possible 
QPO might be present at ~ 0.035 Hz (Fig.[4]i, even if the quality 
of the data does not allow us to firmly establish its presence. On 
the other hand, we cannot exclude that the 0.078 Hz feature is it- 
self a QPO at a frequency close to that of the coherent pulsation. 
Given the faintness of the source, we cannot conclude further on 
that matter. 

In fact, IGR J19294+1816 lies in a region populated by 
Be-HMXBs in the so-c alled "Corbet diagram" dCorbetl 119861: 
ICorbet & Krimml l2009t) , as demonstrated in Fig. [6] Note that 
this plot is the most recent update of the Corbet diagram for 
INTEGRAL sou rces as of June, 2009 ( Bodaghee et al. 2009 in 
prep.). Recently, [Bodaghee et al. (2007) further explored the pa- 
rameter spaces occupied by high-energy sources. They noticed 
that Be-HMXB and supergiant HMXB also segregate in differ- 
ent parts of the Nh vs. orbital period and Nh vs. spin period 
diagrams. The value of the absorption we obtained through our 
spectral analysis, combined with the values of the spin and or- 
bital periods, make IGR J 19294+ 18 16 lie in a region populated 
by Be-HMXBs in these diagrams as well. At first order, one can 
easily understand the segregation in these three different dia- 
grams as results of the age, type of accretion, and probable type 
of orbit. Be systems are younger, and hence have eccentric or- 
bits, longer orbital periods, and shorter spin periods, whereas su- 
pergiant systems are older, are mostly circularised with shorter 



orbital periods, and longer spin periods. In these latter systems, 
in addition, the compact object is embedded in the wind of the 
companion (the feed matter for accretion) which explains their 
(usually) higher intrinsic absorption. In this respect, all our re- 
sults point towards IGR J 19294+ 18 16 being a Be-HMXB (the 
same conclusion, although only based on the Corbet diagram, 
was presented bv lCorbet & Krimml 120091) . 

INTEGRAL has unveiled a new (sub-)type of supergiant 
HMXB, characterised by short (-hours) and intense flares seen 
at X-ray energies: the so-called SFXTs. Models of SFXTs in- 
volve stochastic accretion o f clumps from the (heterogeneous) 
wind of the supergiant (e.g. lin't Zandl 120051 : iNegueruela et all 
20061). on top of longer, but fainter, periods of activity (e.g. 



Sidoli et ail 120071) . In this respect, the detection of periods of 



short and intense flares in IGR J 19294+ 18 16 with INTEGRAL 
(Fig. [B, a behaviour typical of SFXTs, raises the possibility 
that this source belongs to this new class. Its position in the 
Corbet and the Nh vs. spin or orbital period diagrams (Fig. |6| 
is not a definitive contradiction to this hypothesis, since system- 
atic X-ray studies of SFXTs have shown that, contrary to other 
HMXBs, they seem to populate any part of the diagrams. In par- 
ticular, IGR J18483-0311 and IGR Jl 1215-5952, the only two 
SFXTs for which both orbital and pulse periods are known, lie 
in the region of Be-systems in these representations (Fig. [6]). The 
former source has a pulse period of 21.05 s, an orbital period 
of 18 .5 d, and has a B0. 5 la companion (see iRahoui & Chatvl 
I2008L and references therein). The latter has a pu lse period of 
about 187 s, an orbital p eriod of about 165 d (e.g. ISidoli et all 
~* 20091) . and is associated with a B supergiant 



l2007UDucci et al 



(Negueruela et al 



20051) . To reconcile the fact that they are long 



period systems with a rather low absorption, an eccentric or- 
bit is invoked. IGR J 19294+ 18 16 could be the third member of 
SFXTs (with a known pulse period) having a long and eccentric 
orbit. In this case, this may point towards the existence of an evo- 
lutionary link between Be-HMXBs and eccentric SFXTs (Liu et 
al. 2009 in prep.). The quality of the data is not good enough to 
permit us to make any firm conclusions concerning the nature of 
the system, and only an identification of the optical counterpart's 
spectral type will resolve this issue. 
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